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Abstract. Proposed mechanisms for the formation of km-sized sohd plan- 
etesimals face long-standing difficuhies. Robust sticking mechanisms that would 
produce planetesimals by coagulation alone remain elusive. The gravitational 
collapse of smaller solids into planetesimals is opposed by stirring from turbulent 
gas. This proceeding describes recent works showing that "particle feedback," 
the back-reaction of drag forces on the gas in protoplanetary disks, promotes 
particle clumping as seeds for gravitational collapse. The idealized streaming 
instability demonstrates the basic ability of feedback to generate particle over- 
densities. More detailed numerical simulations show that the particle overdensi- 
ties produced in turbulent flows trigger gravitational collapse to planetesimals. 
We discuss surprising aspects of this work, including the large (super-Ceres) 
mass of the collapsing bound cluster, and the finding that MHD turbulence aids 
gravitational collapse. 



1. Introduction 

Coagulation i s the dominant mechanism for the growth of dust grains via van der 
Waals forces ( Dominik &: Tielend 1997) and for the growth of solid protoplanets 



by gravitational binding ( Goldreich et al. 20041 ) . But sticking is difficult in the 



mm-km size range, a s confirmed by an extensive body of experimental work 
( Wurm &: BlumI |2006| ) . Moreover, incremental growth of planetesimals leads to 



the rapid (~ 10^ yr) inspiral of particles near a meter in size. 

An alternative hypothesis ( Safronovl Il969l : [Goldreich &: Wardlll973l . here- 



after GW) proposes that km-sized planetesimals formed from the gravitational 
collapse of smaller solids. This mechanism overcomes the sticking and radial 
drift obstacles in one fell swoop. However stirring by turbulent gas opposes grav- 
itational collapse. GW noted that their dense particle midplar ie would trigger 
a turb ulent boundary layer via Kelvin-Helmholtz instabilities. IWeidenschilling 



( 1980l ) argued that this particle-driven turbulence (an example of drag force 



feedback) would stir up the particle midplane enough to prevent gravitational 
collapse. This appeared to be a fatal flaw of the GI hypothesis, since particle 
settl ing was a sel f-limi ting process. 



ISekival (jl99^ ) and lYoudin &: Shul ( 20021 ) showed that particles could actually 



help their cause of becoming planetesimals. When the surface density of solids 
relative to gas is larger (by factors of few) than solar abundances, then vertical 
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shear is no longer strong enough to overcome the anti-buoyancy of the dense 
midplane layer. The ability of particles to stir themselves is limited. 

However these works did not model the back reaction of drag forces on the 
gas in detail. Indeed most analyses of mid plane Kelvin-Helmholtz instabilities 
(with the exception of Johansen et al.ir2006i ) reduce the particle and gas dynam- 
ics to a simplified set of equations which omit relative motion and prevent a 
detailed examination of the effects of drag forces. 

Section [2?] describes the streaming instability ( Youdin &: GoodmanI l2005l . 
hereafter YG) which takes the opposite approach of neglecting stratification 
and vertical shear to consider two-way drag forces in a self-consistent, if in- 
complete, model. This idealized system shows that particle fe edback triggers 



sponta neous particle clumping. The detailed 3D simulations of iJohansen et al 



( 2007al . hereafter JOMKHY) include vertical stratification, self-gravity, MHD 
turbulence and multiple particle sizes, see fJST] These models show that clumps 
of 15 - 60 cm boulders, augmented by feedback effects, collapse gravitationally 
into bound clusters, which should continue to contract into planetesimals. 



2. The Streaming Instability 



2.1. Linear Growth 

The streaming instability arises in a simplified model of midplane layers in a 
protoplanetary disk (YG). The model analyzes the local, axisymmetric Keplerian 
dynamics of gas and a single particle species in the absence of stratification. The 
gas and solid components are coupled by a linear drag acceleration which satisfies 
Newton's third law by including the back-reaction of drag forces on the gas. The 
only forcing is a constant radial acceleration, representing the global pressure 
gradient, which induces sub-Keplerian rotation of the gas. 

Due to the neglect of stratification and particle settling, the model has a 
well-defined equilibrium state. The sub-Keplerian headwind robs the particles of 
angular momentum and causes an inward (outward) radial drift for solids (gas) . 
The coupled drift velocities are ( Nakagawa et al.lll986l : lYoudin Sz Johansenll2007l . 
hereafter YJ), 
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Uz = Wz = 0, where Ui {wi) are the velocity components of gas (particles) relative 
to a Keplerian orbit, and {x,y,z) are the locally Cartesian radial, azimuthal, 
and vertical coordinates. The velocity scale, tjvk ~ 20 - 50 m/s is set by the 
global radial pressure gradient (see YG or YJ). The system is characterized 
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by two dimensionless parameters: Tg = J7tstop) the aerodynamic stopping time 
tstop normalized to the Keplerian orbital frequency, fi, is a proxy for particle 
size (bigger particles give larger values of Tg); and // = Pp/pg, the equilibrium 
particle-to-gas density ratioQ 

YG found that equilibrium drift solutions are robustly unstable, with linear 
growth for any values of Tg and p. The instability produces particle overdensities 
as a mechanism to communicate the back-reaction of the drag force on the gas. 
Gas density perturbations exist, but with much smaller amplitudes since motions 
are subsonic. 

The robustness of the instability relies on the free energy contained in the 
relative drift of particles and gas. Unlike the two stream instability in plasma 
physics, which is mediated by electric fields instead of dissipative drag forces, 
rotation is required to drive streaming instabilities. Ultimately the instability 
is driven by the work done by the global pressure gradient on the outward 
drifting gas. Thus t he streaming instability differs from (e.g.) the magneto- 
rotational instability ( Balbus Hawlevlll99ll ). which taps the energy released by 



the outward transport of angular momentum due to orbital (usually Keplerian) 
shear. See YJ§5 for an analysis of disk energetics with particle-gas coupling. 

Though the input physics to the streaming instability is rather simple, the 
resulting behavior is quite complex, partly because the motions contain many 
degrees of freedom. Axisymmetric growth requires a "2.5D" analysis, with all 
three velocity components included for each species plus particle density and gas 
pressure (for a total of 8 components). 

The peak growth rate increases with p as the strength of feedback increases. 
One might expect that growth rates would peak near Tg = 1, where most drag 
related processes (e.g. drift speeds and vortex trapping) are maximized. This is 
indeed the case in the "normal" gas-dominated regime, p < 1. However in the 
particle-dominated regime, > 1 a curious behavior occurs. The peak growth 
rate increases as Tg 0. This behavior is seen in Fig. 3 of YG and Fig. 2 of YJ. 
The rapid growth for Tg 1 with ^ > 1 is only possible because it occurs on 
small length scales. 

Numerical Confirmation The linear growth rates derived by YG were con- 
firmed by numerical simulations in YJ. An Eulerian hydrodynamics code (the 
Pencil Code) was used to model the gas, while Lagrangian "super-particles" 
(because a realistic number of particles is impossible to include) represented 
the solid component. Reproducing the linear growth rates by careful seeding of 
velocity and density eigenmodes was useful for calibration (and algorithm de- 
velopment) of the code. This task also confirmed that the analytic treatment of 
particles as a pressureless fluid in YG was a valid analytical approximation. 

2.2. Non-linear Evolution 

The non-linear evolut ion of the streaming i nstab ility into particle-driven turbu- 
lence was explored in I Johansen &: YoudinI ( 20071 . hereafter JY), using the same 



code as YJ. The qualitative behavior was quite diverse, depending on the param- 
eter values. In simulations with Tg = 1 (roughly 20 cm - m sized boulders), the 



^In YJ and JY, e was used in place of ^ for consistency with previous numerical simulations. 
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linear growth phase was followed by an upward cascade to large clumps which 
survived for tens of orbits where particles collected as in a traffic jam. Smaller 
particles, with Tg = 0.1 did not show a pronounced upward cascade and clumps 
were generated and destroyed in an orbital time. A remarkable surprise (in the 
Ts = 0.1 runs) was that the densest particle clumps drifted inwards faster than 
the particles in voids, a reversal of the density dependence for laminar flows. 

The most interesting statistic (for the purposes of planetesimal formation) 
is the degree of particle clumping. The large particle (ts = 1) runs showed 
the greatest clumping, with peak particle overdensities of several hundred. The 
small particle (ts = 0.1) runs produced particle densities up to factors of tens, 
as long as the mass ratio fi > 1. Though the linear instability is very powerful in 
this regime, the shorter clump lifetime produces more modest (but still highly 
non-linear) overdensities. Strongly non-linear clumping is not universal. In the 
small particle, gas-dominated case (/U = 0.2, Tg = 0.1) then particle density 
fluctuations were order unity at most, consistent with the slow growth rates in 
this regime. 

While peak densities are physically relevant (e.g. for gravitational collapse), 
they are always numerically underestimated, as higher resolution and larger 
simulation boxes will find rarer fluctuations. JY also analyzed the transfer of 
angular momentum (which was inward!) and the diffusion of particles. 

3. The Kitchen Sink Simulations 

While the streaming instability provides an idealized demonstration of the rele- 
vance of particle feedback, the simulations in JOMKHY included a much wider 
range of physical processes, hence the tongue-in-cheek label "kitchen sink." 

3.1. Feedback and MHD turbulence 

The first set of simulations in JOMKHY (see their Fig. 1) studied the effect 
of feedback on particle clumping in magneto-rotational turbulence. Simulations 
were fully 3D with ideal MHD. Vertical gravity was included to allow particle 
settling (but was not necessary, or included, for the gas since only the midplane 
region was considered). Without feedback, particles (with uniform Ts = 1) 
reached a peak density of 10 - 30 times the gas density due to a balance between 
vertical settling, pressure trapping, and (on the other hand) turbulent diffusion. 
This concentration factor is > 10'^ relative to uniform mixing with the gas (since 
the mass fraction of solids-to-gas was fixed at the solar abundance of l%)o 

When the back-reaction of the drag force on the gas is included, peak par- 
ticle densities were amplified by a further factor of 10 (to > 100 times the gas 
density). The ability of particle feedback to enhance particle clumping is robust, 
and not just a consequence of the most idealized case of the streaming instability. 
These simulations with feedback included all viable sources of turbulence; ex- 
ternally imposed (here magneto-rotational), and particle driven both by vertical 
shear and streaming instabilities; and still particle clumping is strong. 



^Recall the comment in i|2.2.l that the peak particle densities are underestimates and are most 
meaningful for comparing simulations at the same resolution. 
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3.2. Gravitational Collapse Simulation 

The particle clumping discussed above should readily trigger gravitational col- 
lapse. Particle densities exceeding the gas density by ^ 100 also exceed the 
Roche density, since PRoche/Pg — /{Gpg) ~ 20 (for double the minimum solar 
nebula as in the simulations). To test this expectation, self-gravity was added 
to the model. For greater realism, these runs included a range of particle sizes 
from 15 - 60 cm (ts = 0.25 - 1.0) and a prescription for cooling by inelastic 
collisions was added. 

To lessen the influence of initial conditions, the simulation was run in three 
phases. First, the magneto-rotational turbulence developed for 10 orbits. Sec- 
ond, particles were added with vertical gravity and two-way drag forces (similar 
to the back-reaction simulations in but with higher resolution and a range 

of particle sizes). After 10 orbits, vertical settling, turbulent clumping and 
streaming instabilities produced peak particle densities almost 100 times the 
gas density (slightly less than in ^3.1.1 due to the range of particle sizes), see 
JOMKHY Fig. 3. 

Finally in the third phase, self-gravity and collisional cooling were turned on 
and the system ran for seven more orbital periods (see JOMKHY Fig. 2). The 
densest clump of particles suddenly became a bound cluster of nearly a Ceres 
mass (equivalent to a 500 km planetesimal) . The cluster continued to accrete 
particles at > (l/3)Mc eres psr Orbital period, for a "final" (but still increasing) 
mass of 3.5Mceres after only seven orbits. 



3.3. Surprises 

There are many surprising features of the JOMKHY simulations which go against 
conventional wisdom. The Ceres mass scale of the bound clusters is striking, 
since it exceeds by many orders of magnitude the traditional estimate of a "km- 
sized planetesimal." This mass need not go into a single planetesimal. By anal- 
ogy with Jeans collapse, one might expect further fragmentation due to colli- 
sional cooling during the collapse. Only the initial stages of collapse have been 
modeled so far. 

Also the Ceres mass scale is not a universal result. The mass scale is 
introduced not by self-gravity acting on a smooth background as in GW, but 
rather by particle clumping in the (imposed and particle-generated) turbulent 
flow. Thus the cluster mass depends on particle sizes and the nature of the 

imposed turbulence. 

Indeed the supplementary information to JOMKHY ( Johansen et al.ll2007bl ) 



considers the case of low ionization disks in §12.1, where magetic fields are 
ignored, and only particle-driven turbulence exists. The bound clusters in these 
simulations have masses and clump accretion rates lower by more than two 
orders of magnitude (see Tables 5 and 6 in the supplement). 

This was perhaps the biggest surprise of the work, that imposed turbulence 
abets rather than hinders the gravitational collapse of solids. The non-obvious 
explanation is that turbulent particle concentration has a "snowball effect." 
Particle clumping produced by magneto-rotational turbulence provides seeds 
where particle feedback produces denser clumps than the streaming instability 
(plus vertical shear) alone can trigger. 
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4. Discussion 

Tlie long-standing mystery of planetesimal formation suddenly appears much 
less daunting, if no less fascinating. For decades, the obstacles seemed insur- 
mountable: low sticking efficiencies, rapid radial migration and the inevitability 
of turbulent stirring. Now it is clear that turbulent stirring does not necessarily 
prevent gravitational collapse of solid particles into planetesimals. As discussed 
in ^3.3.\ more turbulence can promote gravitational collapse in some cases! 

The breakthoughs in particle-gas dynamics are related to the role of particle 
feedback. Once thought only to hinder particle settling by triggering Kelvin- 
Helmholz instabilities, we now realize that feedback triggers and augments par- 
ticle clumping. The link between feedback and clumping has many pillars of 
support: the idealized streaming i nstability (YG, YJ, J Y), Kelvin-Helmholz in- 
stabilities with uniform rotation (Johansen et al.ll2006l ). and now 3D stratified 
models of Keplerian disk midplanes (JOMKHY and its supplement). 

Much work and many uncertainties remain. The large initial sizes assumed 
by JOMKHY might be unrealistic, either because coagulation stalls at smaller 
sizes, or if less violent gravita tional collapse occurs first (see work on dissipa- 



primitive undifferentiated meteorites betray no structures larger than mm-sized 
chondrules suggest this concern may be valid. However, particles which are small 
and tightly coupled to the gas disk are less effected by the dynamical clumping 
mechanisms discussed here, at least at the resolutions studied to date. Perhaps 
more surprises are needed. 
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